Caldeira-Leggett Model, Landau Damping, and the Vlasov-Poisson 

System 1 



George I. Hagstrom and P. J. Morrison 

Department of Physics and Institute for Fusion Studies, The University of Texas at Austin, Austin, TX 

78712. 



Abstract 

The Caldeira-Leggett Hamiltonian (Eq. flTJ below) describes the interaction of a discrete 
harmonic oscillator with a continuous bath of harmonic oscillators. This system is a stan- 
dard model of dissipation in macroscopic low temperature physics, and has applications to 
superconductors, quantum computing, and macroscopic quantum tunneling. The similar- 
ities between the Caldeira-Leggett model and the linearized Vlasov-Poisson equation are 
analyzed, and it is shown that the damping in the Caldeira-Leggett model is analogous to 
that of Landau damping in plasmas [1]. An invertible linear transformation 0, 0] is pre- 
sented that converts solutions of the Caldeira-Leggett model into solutions of the linearized 
Vlasov-Poisson system. 

Keywords: Landau damping, continuum damping, quantum dissipation, Lie-Poisson 
bracket, Caldeira-Leggett, Hamiltonian 



1. Introduction 

In 1946 Landau [l| theoretically predicted the collisionless damping of the electric field 
in a plasma governed by the Vlasov-Poisson system. This result has been of great impor- 
tance in the field of plasma physics, and indeed collisionless or continuum damping, as it is 
sometimes called, occurs in a wide variety of kinetic and fluid plasma models that possesses 
a continuous spectrum. For example, such damping occurs in the context of Alfven waves 
in magnetohydrodynamics (see e.g. Chap. 10 of |4|) and has been proposed as a mechanism 
for plasma heating in response to electromagnetic waves. 

Many other systems also undergo Landau damping, both inside and outside of plasma 
physics. It is not surprising that Landau damping exists in stellar dynamics governed by the 
Jeans equation [|§] because this equation is of Vlasov type but with an attractive interaction 
potential. In fact, Landau damping occurs in collisionless kinetic theories with a rather 
large class of potentials, and recently has been proven rigorously to exist in the nonlinear 
case 0, 0] • Landau damping exists in the context of the fluid mechanics of shear flow (see 



1 Dedicated to Steve Childress whose unique combination of mathematical and physical insight has been 
both an inspiration and a guiding light. 

Preprint submitted to Physica D September 1, 2010 



e.g. 0, [i| which contains a list of original sources over a period of more than 50 years) 
and the description of wind driven water waves. It also appears in multiphase media 10 
and has been established for systems containing large numbers of coupled oscillators, most 
notably the Kuramoto model. This has implications for biological models describing the 
synchronization or decoherence of the flashing of fireflies and chirping of crickets as well as 
other phenomenon in mathematical biology 

Another class of continuum systems involves the interaction of a discrete oscillator with 
a continuous bath of oscillators. In these systems the oscillator can be a particle or one 
mode of some field, and the bath often represents thermal fluctuations or radiation. One of 
the first detailed treatments of such a system is due to Dirac [12], but early on Van Kampen 



also used such a model to describe the emission and absorption of light by an atom [13| . The 
single wave model of plas ma p hysics, which describes both beam plasma and laser plasma 
interaction physics 1% [15|, llyL is al so an example. The example of interest in this paper is 



the Caldeira-Leggett model [17 . 

The Caldeira-Leggett model was invented in order to study quan tum tunneling in the 



presence of dissipation and the quantum limit of Brownian motion [18(. A model of this type 
was deemed necessary because quantum mechanics is incompatible with frictional forces. 
However, the Caldeira-Leggett model is a Hamiltonian system that exhibits dissipation by 
coupling to a continuum, i.e., it has Landau damping. The Caldeira-Leggett Hamiltonian is 
the sum of the Hamiltonian of a classical harmonic oscillator, the Hamiltonian of continuous 
bath of harmonic oscillators, and a linear coupling term between the discrete and continuous 
degrees of freedom. The discrete degree of freedom corresponds to a macroscopic system and 
the bath of oscillators represent the environment. The coupling causes the discrete oscillator 
to damp by transference of energy to the continuum. This system has become a standard 
model for studying the physics of low temperature quantum systems, and it has numerous 
applications ranging from the understanding of superconducting circuit elements to qubits 
in quantum computers 
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We analyze the classical Caldeira-Leggett model using a procedure analogous to that 
used by Landau to analyze the Vlasov-Poisson system of plasma physics. Following Landau, 
the initial value problem can be solved using the Laplace transform and the rate of decay 
can be derived in the weak damping limit. This paralleling of Landau's original calculation 
suggests a connection between this system and the Vlasov-Poission system. In fact, we will 
show that both systems can be mapped into a normal form that is common to a large class 
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snow tnat Dotn systems can De mapped into a normal iorm tnat is common to a large 
of infinite-dimensional Hamiltonian systems that have a continuous spectrum 0, 0, @, 

The Caldeira-Leggett model, like all Hamiltonian systems in the class, has a continuous 
spectrum that is responsible for the damping through phase mixing (filamentation) and 
the Riemann-Lebesgue lemma. Because this structure is shared by a number of important 
physical systems, it is interesting to determine the nature of their similarities. It is well- 
known that the properties of linear ordinary differential equations are closely tied to the 
spectra of their time evolution operators. In fact, for given spectra there are a number 
of normal forms. Any linear finite-dimensional Hamiltonian system can be reduced to one 
of these normal forms (ODEs) through an appropriate transformation, and in this sense 
the behavior of such systems is completely understood. The theory of normal forms for 
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infinite-dimensional Hamiltonian systems is not nearly as well-developed as that for finite- 
dimensional systems, but for some systems much is known. For systems with continuous 
spectra, the analog of diagonalization is conversion into a multiplication operator. If the 
original system is / = Cf, then a transformation T such that TCT' 1 is a multiplication 
operator would diagonalize the system. Any two systems that have the same normal form 
would thus be equivalent through some linear transformation. 

This procedure has been performed for the linearized Vlasov-Poisson equation jjjj, and 
when the spectrum is purely continuous the time evolution operator is equivalent to the 
multiplication operator x. This discovery led to the discovery of an entire class of transfor- 
mations diagonalizing linear infinite-dimensional Hamiltonian systems of a certain form 20 . 



In fact, it is always possible to perform such a transform in the special case of a bounded, self- 



adjoint operator [21] . The operators dealt with here are usually unbounded and non- normal 
(even if they did exist in a Hilbert space), as is often the case when dealing with contin- 
uous Hamiltonian matter models. A precursor to the discovery of such transformations is 
existence of a complete basis of singular eigenfunctions of the original equation, a treatment 



that is common for systems with continuous spectra that dates back to Dirac 12] . In fact 



these methods have been developed in parallel within the field of plasma physics beginning 
with the work of Van Kampen 22] and within condensed matter physics through the work 
of Dirac and later Fano [23]. Caldeira and his collaborators developed a diagonalization 
method for the Caldeira-Leggett model j24j], although they were primarily interested in the 
time evolution of the discrete degree of freedom and thus did not write down the full inverse 
of their transformation. In this paper we complete the treatment of the Caldeira-Leggett 
system. Then, we note that the normal form is the same as that of the Vlasov-Poisson 
system and that the models are thus equivalent through the use of an integral transform. 

Specifically, in Sec. [2] we review the Caldeira-Leggett model and then, in the spirit of 
Landau PJ, present its Laplace transform solution in Sec. 21 This is followed by obtaining the 
singular eigenfunctions, in the spirit of Van Kampen [22] and Dirac 12], and the invertible 
integral transform akin to that of [3] for transforming to normal form. In Sec. [5] we show 
explicitly how the Caldeira-Leggett model is equivalent to a case of the linearized Vlasov- 
Poisson system. Finally, in Sec. |6] we conclude. 



2. Caldeira-Leggett Hamiltonian 



As noted above, the Caldeira-Leggett model is an infinite-dimensional Hamiltonian sys- 
tem describing the interaction of a discrete degree of freedom with an infinite continuum of 
modes 18] . The continuum is typically referred to as the environment. The Caldeira-Leggett 
model has the following Hamiltonian: 



HcL[q, p; Q, P] = §P 2 + \ (n + £ dx Q 2 



+ dx 



x 



^(p(x) 2 + q{xf) + Qq(x)f(x) 



(1) 
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which together with the Poisson bracket 



fdAdB dAdB\ f (5A5B 5A5B\ 

1,1 ~ [dQdP ~ dPdQj + J^ x ~ ^) (2) 

produces the equation of motion for observables in the form F = {F, H}, where F is any 
functional of the discrete, (Q,P), and continuum, (q,p), coordinates and momenta. Note, 
it is assumed that f(x) is chosen so that the integrals of ([1]) exist. The coefficient of Q 2 
includes a frequency shift term that is used to make the Hamiltionian positive definite. We 
take p and q to be functions on the positive real line, R+, and P and Q to be real numbers. 
Hamilton's equations for the Caldeira-Leggett system are thus, 

q(x) = xp(x) (3) 
p(x) = —xq(x) — Qf(x) (4) 
Q = QP (5) 

P = - (tt + [ dx Q- I dx q{x)f{x) . (6) 

V ./to , 2x I ./» 



This system was originally introduced by Caldeira and Leggett in 1981 [12J. They initially 
considered a very massive harmonic oscillator coupled to a large number of light harmonic 
oscillators with varying frequencies, and then studied the limit of the light oscillators becom- 
ing a continuous spectrum. The coupling causes Q to decay to zero with time, and therefore 
the system can be used to model dissipation. This makes it an ideal system to model the 
effects of dissipation in quantum mechanics and especially quantum tunneling. It has been 
extensively studied and is frequently mentioned in the condensed matter literature. There 
have been some controversies about the physics of the damping and the physicality of the 



initial conditions [19j. Connecting this system with plasma physics, where much intuition 
has been developed over the years about wave-particle interaction, can help to improve the 
understanding of its behavior. For example, a clear picture of filamentation can be viewed 
in the numerical work of (25l |. 

Systems with continuous spectra exhibit phase mixing or filamentation. A wide variety 
of systems have this property, and of course most famously the Vlasov-Poisson system for 
which Landau first discovered his damping. We will demonstrate explicitly that the damping 
mechanism of the Caldeira-Leggett model is Landau damping. Furthermore we will derive 
a transformation that converts the Caldeira-Leggett model into the equation describing the 
time evolution of a single mode of the linearized Vlasov-Poisson system. Two copies of 
the Caldeira-Leggett model will be canonically equivalent to the time evolution of a pair of 
opposite k Fourier modes of the linearized Vlasov-Poisson equation. 

3. The Landau Solution and Landau Damping 

One of the classical calculations in plasma physics is the solution of the linearized Vlasov- 
Poisson equation using the Laplace transform. This yields a formula for the solution of the 
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initial value problem and also facilitates the derivation of the damping rate for the electric 
field. It is possible to do the same thing for the Caldeira-Leggett model. We begin with the 
set of Hamilton's equations that were written down in the previous section and eliminate 
the two momenta to derive a pair of second order equations for the coordinates, 



q(x) = —x 2 q(x) — Qxf(x) 

Q = -n 2 c Q-tt [ dxf(x)q(x), 



where for convenience we use the corrected frequency, 

tt 2 c := tt 2 + tt [ dxf(x) 2 /2x. 
Jr + 

Defining the Laplace transform of the coordinates by 



q(x, s) = dt q(x, t)e 



St 



Q(s) = / dtQ(t)e 



-St 



results in the following set of algebraic equations: 

s 2 q(x, s) = —x 2 q(x, s) — Q(s)xf(x) + sq(0, x) + g(0, x) 
s 2 Q(s) = -tt 2 c Q(s) -SI [ dxq(x, s)f(x) + sQ(0) + Q(0) , 



which can be easily solved for Q(s) 



Q(s) = 



■J 



ft I dx 



f(x)(sq(0,x) + q(0,x)) 



S 2 + X 2 



s z + n z -n / dx 



+ sQ(0) + Q(0) 
xf(x) 2 



S 2 + X 2 



The Laplace transform is inverted using the Mellin inversion formula, 

Q(t) = — / dsQ(s)e st , 

2,1X1 Ja_ ir ~ 



(7) 
(8) 



(9) 



(10) 
(11) 

(12) 
(13) 



(14) 



(15) 



where (3 is any real number that ensures Q(s) is analytic for Re(s) > (3. This integral is 
usually evaluted using Cauchy's integral formula, whence asymptotically in the long-time 
limit the behavior of the solution is given by the poles of Q(s). Thus, the solution will be 
dominated by an exponentially decaying term arising from the pole of Q(s) closest to the 
real axis. We assume the closest pole is indeed close to the real axis and that there are no 



5 



poles with a positive real part, i.e., that the solutions are stable. This is the weak damping 
limit. As long as f(x) is Holder continuous, the poles of Q(s) come from the zeros of the 
denominator, so we are interested in the roots of the equation 



Q = s 2 + n 2 n f dx ^L = s 2 + n 2 n f dx {(M)- , (16) 

Jr + s 2 + x 2 J r 2{x-%s) 

Here /(Ix)) 2 . is the antisymmetric extension of f(x 2 ) defined by = sgn(x)f(\x\) 2 . 

Making the substitution uj = is, yields the dispersion relation, which in the limit uj tends to 
the real axis becomes 

2 n^ Q Ij /(H)- i M t(\ h2 n tn\ 
u -Q c + -fdx + _/(ar)_ = 0, (17) 

2J R x-u 2 

where j- denotes the Cauchy principal value integral. For quantities not yet integrated, we 
will denote this by PV. This equation can be viewed as the dispersion relation in the weak 
damping limit. Let u c be a real solution to the real part of the above equation. Then let 
u be a root of the previous equation, assume 7 = Im(u) is small, and solve for 7 to first 
order; i.e. = 2iu) c ^ + inflf(\uj c \) 2 _/2 or 

7 = -^V(M) 2 - (18) 

There are a large number of methods used to derive damping of Q for this model. The 
standard approach is to attempt to prove that after suitable approximations Q satisfies the 
equation of motion of a damped harmonic oscillator. The treatment here is almost identical 
to the method that was used to treat the Vlasov-Poisson equation by Landau, and agrees 
with other derivations of the damping rate in the weak damping limit [17|. 



4. Van Kampen modes: diagonalization of the Caldeira-Leggett model 

The Laplace transform method is just one way to treat the Vlasov equation. Another 
way is to write the solution as a superposition of a continuous spectrum of normal modes, a 
method attributed to Van Kampen [22j. Such modes of the Vlasov equation are called the 
Van Kampen modes, and we will see that they exist for the Caldeira-Leggett model as well. 
We formally calculate the Van Kampen modes for this system and use them to motivate 
the definition of an invertible integral transform, akin to those of 0, 0, 0, |2(J, that maps 
the Caldeira-Leggett model to action-angle variables, the normal form for this Hamiltonian 
model. The nature of the transformation depends on the coupling function f(x). In the 
present treatment we will assume that /(0) = 0, but that / does not vanish otherwise. 
We will also assume that the dispersion relation does not vanish anywhere. This excludes 
the possibility of discrete modes embedded in the continuous spectrum. The case where 
the Caldeira-Leggett model possesses such modes will be treated in future work. As stated 
above, the normal form of the Caldeira-Leggett Hamiltonian will be seen to be equivalent to 
that for the Vlasov-Poisson system through the integral transformation introduced in 0, IsJ. 
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The first step is to obtain a solution with time dependence exp(— iut) and derive equations 
for the amplitudes of a single mode (q u ,Pu, Qu, P u )- To this end consider 

iuq u (x) = -xp u (x) 
iup u (x) = xq u (x) + Q u f{x) 

iuP u = (n+ [ dx Qu + [ dxq u (x)f(x). (19) 

V Jr+ ^ x J Jr+ 

Note, although we use the subscript, u G R is a continuum label. Eliminating the momenta 
from Eqs. (jl~9]) yields 

(u 2 - x 2 )q u (x) = Q u xf(x) (20) 
(u 2 - n 2 c )Q u = Q [ dxq u (x)f(x) , (21) 



where recall Q c is defined by Q. Of these, f[2"0l is solved following Van Kampen (a gener- 
alized function solution that dates to Dirac |12J) giving the general form for q u 

q u {x) = PV 9flM + C u Q u S(\u\ - x) . (22) 
u z — X A 

Substitution of (1221) into (12T1) determines C u , 



u 2 -n 2 c = ni dx + nc u f(\u\) (23) 



u 2 — X 2 

,2 _ n2 /• f(\„\\1 



C u2 - n < Lx m) ~ (0 1) 

u n/(|u|) j/ x 2( M -x)/(| M |)- l ' 4j 

Therefore we can specify an initial condition on the amplitudes Q u and compute the cor- 
responding coordinates and momenta by an integral over the real line. Each mode oscillates 
with a different real frequency, with the expression for the solution given by 

q(x,t) = ldu Q f f{ - X } e~ m + / duC u QJ(\u\ - x)e~ mt (25) 

Q(t)= \duQ u e~ lu \ (26) 

were Q u acts an amplitude function that determines which Van Kampen modes are excited. 

The Caldeira-Leggett model can be diagonalized and solved by making use of the integral 
transform alluded to above. Previously, Caldeira et al. [24| derived a transformation that 
diagonalizes the Caldeira-Leggett model. However, they were interested in solving for the 
evolution of the variable Q and therefore did not attempt to write down the full inverse of 
the operator (except in a special case where they made use of the evolution of the reservoir). 
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We will extend their results by deriving the inverse map that we will use to establish the 
equivalence with the Vlasov-Poisson system. 

In order to define the transform, we introduce a number of other important maps and 
introduce our notation. Extensive use will be made of the Hilbert transform, which is defined 
for a function g(x) on R by 



»[«](«) 



1 

* Jm 



dx 



x — V 



We also need some Hilbert transform identities 26, 0]. Let g, g±, and g2 be functions of 
iGt and suppose that all the expressions we write down are well defined, then the following 
hold: 



H[H[g]] = -g 
H[ 9l H[g 2 ] + g 2 H[ gi }} = H[ 9l ]H[g 2 ] - g x g 2 

H[vg] = vH[g] + - 

7T 

Next we define two functions, €r and ej by 



dx g . 



6l = irf(xf 



and 



x 



it 



+ nH[f(\x\ 



(27) 
(28) 

(29) 



(30) 



These together with \e\ 2 := ej + e 2 R are used to define the following integral transforms: 
Definition For functions h(x) on R + , the transform 

T + [h](u) := e R h(\u\) + erH[h(\x\)](u) , 

while 

T + [h](u) := ^h(u) - -^H[h(\x\)](u) , 
\ e \ \ e \ 

Related to the above transforms are two more transforms, 

T.[h]{u) := e R h(\u\) + ejHisgnWhdxD^u) , 

and 

T_[h](u) := ^h(u) - ^H[sgn(x)h(\x\)](u) . 
\ e \ \ e \ 

Using the transform T + it is possible to write the map from the amplitudes of the Van 
Kampen modes Q u to the functions (q(x), Q). To see this consider the expression for q(x) in 
terms of the amplitude function Q u , and simplify it using the Hilbert transform as follows: 



q{x) 



du Q^1Sl!L- -\- I du C u Q u 5(\u\ — x) 



du 



u 2 — x 2 
xf{x)Q u 



2u 

TTxf(x) I H 





1 

i 


1 


1 1 

\u — X u 


+ X 


Qu 

2u 


(x) + H 


Qu 

2u 



C^xQx C_ X Q_ X 
+ 2C X (Q X + Q„ 



-x 



(31) 



Next, decompose Q u into its symmetric and antisymmetric parts: Q u = Q +u + Q- u and 
observe that the antisymmetric parts vanish from both sides of (I3ip . 

q(x) = irxf(x)H[Q +u /v](x) + 2C X Q +X 

= nf(x)H[Q +u ] + (2 - ldx> f{lX Y- ) Q +x , (32) 

where the second line follows from the third Hilbert transform identity combined with the 
fact that Q u+ /u is antisymmetric and thus has a vanishing integral. Now multiply both 
sides of f[3"2l by f(x) and find 

f(x)q(x) = nf(x) 2 H[Q +u ] + (2 - jdx' j^^j Q +x 

= £iH[Q +u ] + e R Q +x 

= T+[Q +U }. (33) 
Now we are set to define a transformation. 
Definition Let Q +u be a function on R + , then the map 



h[Q+] ■= (-^-T + [Q u+ },2 [ du Q uA 
\J[ X ) Jm+ 



The map I C [Q + ], a map from the Van Kampen mode amplitudes to the original dynamical 
variables, has an inverse. To see this note that €r = S + H[ei], where S = 2(x 2 — fll)/fl, 
and let g be a function on M + . Then, 



H[Sg(\x\)] = SH[g{\x\)) + — / dx g , 

Jr+ 



where we have used our Hilbert transform identities to move the x 2 outside of the Hilbert 
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transform of g. Using this, consider the following sequence of identities: 

t + [t + M. 



^(e R g + ejH[g]) - ^H[e R g + e 7 tf [<?]] 



~R 



g + e -^H[g] - ^H[g] - ^[H^g + ejH[g]] - 



dx g 



%9 + - e -^H[g] - ^(HlejWg] - ge,) - / dx g 



e R e I 



H[ei]H\g] 



€1 Au 



dx g 



6i Au 



dx g 



= 9 



6j Au 



dx g , 



(34) 



where in each step use has been made of the various identities above. Because the integral 
of Q + is equal to Q/2, we can define the inverse of / as follows: 

T c [ q ( x ),Q]=Tl[f(x)q(x)] + ^^Q. 

The above transform ignores the (p, P) variables and only produces the symmetric part 
of the Van Kampen modes. We derive the other half of the transformation from a mixed 
variable generating functional. To this end, define Q + = Q and Q_ = P and rescale the 
coordinate part of the transformation by choosing Q = 2 J R+ du Q^6i / {ix\e\ 2 ) : 

T[q(x),Q] 

= i[Q] ■ 



Q 

(Q,q(x)) 



7T 



T + [f(x) q (x)] + ^^-Q 



7r e 



Q 



(35) 
(36) 



Then we introduce the mixed variable type-2 generating functional 



duPI[q(x),Q], 



which produces the transformations in the usual way: 



p{x) = 5 -^ = f{x)T + ] 



7r e 



-P 



pJ-L 

5Q 



du 



2uP 
IT 



7r e 



(37) 
(38) 
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Calculating the adjoint of T simplifies the resulting expression for p(x), viz., 

1 „ 



p(x) 



7r e 



(39) 



Now, analogous to / we define the operator J[P] = (p(x),P) by 

2uP rrr 



j[p\ 



^p 

_l 12 

7r e 



V 7r|e 



(40) 



which can be inverted through the use of the Hilbert transform identities. Define P c = 
P^/n\t\ 2 / 1[, and consider the expression 

J^P( X ) ~ J^2 H [ s 9n(x)p(\x\)] 

= ^(eiH[sgn(u)P c ] + e R P c ) - ■^rH[sgn(x)e I H[sgn(u)P c ] + sgn(x)e R P c ] . (41) 



Paralleling the method used to invert the map from Q to (q,Q), we see a difference 
occurs when evaluating the term ejH[sgn(x)eRp]/\e\ 2 , i.e. 



du uP r 



-^H[sgn(u)u 2 P c ] = x 2 -^H[sgn(u)P c ] + -t-^ f l( „ 

I e l l e l 71 l e l 

= x^ 2 H[sgn{x)P c ] + -^-P. 



(42) 



With this expression we can directly use the inversion calculation for the (q, Q) case to 
obtain the following expression for the full transformation: 



P = 



Q = 



7r e 



T4f(x)p(x)} + 



2 e 7 



7r e 



7T 



r;[/(a;)9(x)] + ) . 



(43) 
(44) 



Applying this transformation to Hamilton's equations yields the equations for a contin- 
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uum of harmonic oscillators. This can be seen directly for both P and Q. 



Q 



7r e 



7T 



2u 61 



T+[xf(x)p(x)] 



P 



it e 



7T 



7r e 



uf(u)p(u) 



uf(u)p(u) 



^H[\x\f(\x\)p(\x\)] + ^^P 



^uH[sgn(x)f(\x\)p(\x\)] + -^-P 
er vr er 



= hP. 
Similarly, for P, 



(45) 



7r e 



7r e 



r:[/(x)p(x)] + 



7T 



7T e 



7T 



ey 



2 6! 



-uT + [f(x)q(x)] + -j^ I dxf(x)q(x)-T4f(x) 2 ]Q 
2il s ej 



7T 



7T 



ey 

= -uQ . 



-uT + [f{x)q{x)} 
-uT + [f(x)q{x)] 



f(x) 2 Q + -^H[sgn(x)f(\x\) 2 ]Q 



2tt s 6/ 



7T 



Q 



2u 2 ti 



Q 



(46) 



Let this full map be called If and consider If as a map from the Banach space L p xL p xM. 2 , 
p > 1, to the Banach space L p x L p . The operator If is a bounded linear functional between 
these two spaces, because each term is either a multiplication operator that is bounded on 
L p , an L p function, or a bounded function multiplied by the Hilbert transform, which is 
another bounded operator. In order to establish the equivalence with the normal mode it 
is important to specify the phase space of the dynamical variables. Using this map we can 
simply choose each functional space to be L p and have a well defined map in each case. This 
map demonstrates how the Caldeira-Leggett model can be written as a superposition of a 
continuous spectrum of singular eigenmodes. 
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Because the transformation to the normal form was a canonical one, the normal form 
Hamiltonian should be the original Hamiltonian of the Caldeira-Leggett model written in 
the new coordinates. We will verify this by direct substitution. For convenience we introduce 
the quantities 

A = —P 2 + - dxxp(x) 2 
2 2 J R+ 

b = t q2+ f (h {x)2 + f ^ q ^ Q ) dx ' 



where £l s = Q 2 C /Q. Evidently, H CL = A + B. Then, 
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Now, analyzing the entire expression for B in a sequence of steps 

n 



B 



dx Qf(x)q(x) 



du Qm 



'7r e 



e_ 



7T e 



T + [/(x)g(x)] + yQ 2 



duQu\ T + [f(x)q(x)] + Q s Q \ duQ. + — / dxT. 



7T e 



7r e 



■0 



+ 



7r e 



-7+ [/(*)?(*)] 



2« 2 = 



Cte ( ( — Q + 7Tif[^(x)/(|x|) 2 ]g 



7r e 



" -,Q] 



7r e 



- j du uQ 2 . 



(49) 



With (jUD and @]) we obtain if CL = J R duu (Q 2 + P 2 ) /2 - the Hamiltonian for a 
continuous spectrum of harmonic oscillators and the normal form for the Caldeira-Leggett 
model. 



5. Equivalence to the Linearized Vlasov-Poisson Equation 



The treatment of the Caldeira-Leggett model of Sec. H] is similar to an analysis of the lin- 
earized Vlasov-Poisson equation performed in 0, [27|. In those papers an integral transform 
was presented that transforms the Vlasov equation into a continuous spectrum of harmonic 
oscillators. The two systems are identical except the spectrum of the Caldeira-Leggett model 
only covers the positive real line. Now we explicitly produce a transformation that takes 
one system into the other. 

The Vlasov equation describes the kinetic theory of a collisionless plasma. Spatially 
homogeneous distribution functions are equilibria, and linearization about such states are 
often studied in plasma physics. In the case of one spatial dimension and an equilibrium 
distribution function fo(v), the linearized Vlasov-Poisson equation around fo is given by 



dt dx 
d 2 ch 



_____ f' 

m dx 







dx 2 



-Aire I dv f 



(50) 
(51) 



where f' Q = dfo/dv. These equations inherit the noncanonical Hamiltonian structure of the 
full Vlasov-Poisson system 28] and have a Poisson bracket given by 



{F,G}j 



dxdv fo 
14 



5F_ 5G 
~Sf'~6f 



(52) 



This bracket is of a form that is typical for Hamiltonian systems describing continuous media 
(cf. e.g. 29, [3^]). The Hamiltonian is given by 

*L = ~I [dvdxv?- + -L [dJ^-X , (53) 



2 J J fo 8tt J \dx 



and the Vlasov-Poisson equation can be written as / = {/, H L } L . 

The spatial dependence of the Vlasov-Poisson equation can be removed by performing a 
Fourier transform. This allows the potential to be explicitly eliminated from the equation 

° lk - i*vj k - ^-f (v) [dvf k = 0. (54) 



dt mk 

The Hamiltonian structure in terms of the Fourier modes has the new bracket 

{F , G]L = ±*[ dv}i (£MMlZ\. (8B) 

^mJi \°fk Sf-k of k 5f-kJ 

and the Hamiltonian functional is simply ( 155]) written in terms of the Fourier modes. 
One way to canonize this bracket is with the following scalings: 



qk(v,t) = f k and Pk(v,t) = ^rr^r 7 (56) 



where k > 0. In terms of these variables the Poisson bracket has canonical form, i.e. 

'5F5G 5G5F' 



{F,G} L = jr [dv( 
fc=i v 



5qk fyk 5q k Sp k 



(57) 



From this point it is possible to derive a canonical transformation that diagonalizes the 
Hamiltonian. We make the following definitions: 

^-^Stk ^W = l + *M. (58) 

G t [f]=E R f + E,H{f] G k [f\ = 5*f-lL.H\f\. (59) 

It was proven in ^ that Gj. = G k ■ 

A transformation to the new set of variables (Qk,Vk) that diagonalizes the system will 
be given in terms of the variables {qk,Pk)- To this end we first introduce the intermediate 
variables (Q^,^) defined by 

qu = G k [Q' k ] and Q' k = G k [q k } . (60) 
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The corresponding momentum portion of the canonical transformation is induced by the 
following mixed variable generating functional: 



oo „ 

H<lk,V' k ]=J2 duP' k G k [q k ], (61) 
k=i Jr 

whence we obtain via Q' k = 8 J- '/SV' k and p k = 8T j8q k , 

£ k = G~ k {q k ] and V' k = G k \p k ] (62) 
Then, the variables (Qk,Vk) are defined as 

Q k = (Q' k -iV' k )/V2 and V k = (V k -i& k )/y/2, (63) 
in terms of which the Vlasov-Poisson Hamiltonian has the form of a continuum of harmonic 



oscillators (see [27| for an explicit calculation), 

oo „ 

H L = J2 du(<£ + Vl)/2. (64) 
k=i J * 

Thus, for a single value of k, this is the same normal form as that of the Caldeira-Leggett 
model, with the exception that the integral here is over the entire real line instead of just the 
half line. If we consider two copies of the Caldeira-Leggett model the normal form would be 
the same as that for a single k value of the linearized Vlasov-Poisson system. By composing 
the transformation that diagonalizes the Caldeira-Leggett model with the inverse of the 
transformation that diagonalizes the Vlasov-Poisson system we obtain a map that converts 
solutions of one system into solutions of the other system. Explicitly suppose that we have 
two copies of the Caldeira-Leggett Hamiltonian, with the same coupling function f(x). Then 
set the normal form of the second copy equal to the normal form of the Vlasov equation on 
the negative real line. Let (qx(x),pi(x), Qi, Pi) be one set of solutions to the Caldeira-Leggett 
model and let (q2(x) , p 2 (x) , Q 2 , P 2 ) be another and let Q(x) be the Heaviside function. Then 
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we can write a solution to the linearized Vlasov-Poisson equation using the following map: 



f k (v,t) = G k 



V2 



Q(u) 



7r e 



2u €i 



/7r|e|(— -u) 2 / -~ rn . , . , — 2-u er(— m) ^ 



€i(-u) 



ttQ \e(-u)\ 2 



+ iG(u) 



+ iG(-u) 



m 



-iG(u) 
-iG(-u) 



'7T e z 



2 e 7 



T_[/(x)p 1 (x)] + - r ^Pi 



7r e 



7T e — u) 



CA-U) 



T4f(x)p 2 (x)} + 



2 eA-u) 



■Pi 



it \e\—u] 



7r e 



2-u 6/ 



7r e — tt 



7r e z 



6/ 



T4f(x)p 1 (x)] + -^-P 1 



tt e 



7T el — M 1 



e/(-u) 



T_[/(x)p 2 (x)] + 



2 €/(-«) 



7T e — M) 



(65) 
(66) 
(67) 
(68) 
(69) 
(70) 
(71) 
(72) 



This map is invertible using the formulas presented earlier in the paper. Given a single 
mode of the linearized Vlasov-Poisson system, fk(v,t), we can write two solutions to the 
Caldeira-Leggett model as follows: 



(q 1 (x,t),Q 1 (t)) = I[-^G[f k ](u,t) +G[f k ](-u,t))} 
( Pl (x,t),P 1 (t)) = J[^(G[f k ](u,t) - G[f k ](-u,t))\ 
Mx, t), Q 2 (t)) = I[^(G[f k ](u, t) + G[f k ](-u, t))] 
(p 2 (x,t),P 2 (t)) = J[^(G[f k ](-u,t) - G[f k ](u,t))\ 



(73) 
(74) 
(75) 
(76) 



Therefore one would expect the solutions of the Caldeira-Leggett model to share the 
same properties as the solutions of the Vlasov-Poisson system. 

It was remarked earlier that both systems exhibit damping. In the Vlasov-Poisson case 
the electric field decays, and in the Caldeira-Leggett model it is the discrete coordinate Q. 
The existence of the transformation between the two systems gives us a way to understand 
what determines the damping rate in each case. In the standard calculation of the Landau 
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damping rate for the Vlasov equation, it is clear that the rate depends only on the location of 
the closest zero in the lower half complex plane of the dispersion relation, which only depends 
on the equilibrium /q. The same is true for the Caldeira-Leggett model, where the damping 
of Q depends on the coupling function /. It is clear that integral transformations change 
the rate of damping, as all the instances of the Vlasov equation and the Caldeira-Leggett 
model share the same normal form but generally have different damping rates. 

It is possible to interpret Landau damping using the normal forms and canonical transfor- 
mation. The dynamical variables of the normal form have a time evolution ~ exp(— iut). The 
observables can then be expressed as some operator on this oscillatory dynamical variable. 
The result will be an oscillatory integral over the real line, and by the Riemann-Lebesgue 
lemma we know that that such an integrated quantity will decay to zero in the long-time 
limit. For the systems at hand, this integral can be deformed into the lower half complex 
plane, and Cauchy's theorem can be used to see that the behavior is governed by the lo- 
cations of the poles of the analytic continuation of the oscillatory integrand. These poles 
determine the exponential damping rate. In these systems the poles are clearly introduced 
by the continuation (following the Landau prescription) of the dispersion relation in the inte- 
gral transformations, which is therefore the origin of Landau damping. We will demonstrate 
this explicitly for the damping of the coordinate Q in the Caldeira-Leggett model. 

Starting from the solution, 

Q(t) = I du (Q(|u|) cos(ut) + sgn{u)P{\u\) sin(wt))^^ 



du (I[q(x), Q] cos(wt) + J[p(x), P}sgn(u) cos(wt))^p , (77) 

we see that each term in the integrand of (1771) has an oscillatory part and has poles at 
the zeros of |e| 2 . The damping rate will be based on the closest zero of |e|, the dispersion 
relation for the Caldeira-Leggett model. Likewise, for the Vlasov-Poisson system we can 
write a similar expression for the density Pk{t), 



p k (t) = dvG k [G k [f]e- lut ] 
Jr 

= [du (e R (G k [f}e- lut ) - HleAGkifle-* 
Jr v 

duG k [f]e- iut . (78 



The damping rates are given by the poles of G k , and the observed rate will be due to the 
closest zero of \e\ 2 to the real axis. Therefore, mathematically the source of the damping 
in the Vlasov-Poisson and Caldeira-Leggett models are identical, it being the nearest pole 
introduced by the integral transformation that diagonalizes the system. 



18 



6. Conclusion 



To summarize, we have shown how the Caldeira-Leggett model can be analyzed the 
same way as the Vlasov-Poisson system. We wrote down the solution using the Laplace 
transform, an expression for the time evolution as an integral in the complex plane over the 
initial conditions. It was then indicated how Cauchy's theorem can be used to derive the 
time asymptotic behavior of the solution, and it was described how the long-time damping 
rate is equal to the distance from the real axis of the closest zero of the dispersion relation 
(when analytically continued into the lower half complex plane). Thus, the damping of the 
Caldeira-Leggett model can be seen to be a rediscovery of Landau (or continuum) damping. 
Caldeira and Leggett introduced their system to study damping in quantum mechanical 
systems, and it is now seen to be one of many interesting physical examples of Hamiltonian 
systems that exhibit such behavior. 

Next we described how to analyze the Caldeira-Leggett model by means of singular 
eigenmodes, paralleling Van Kampen's well-known treatment of the Vlasov-Poisson system. 
Here the solution was written as an integral over a distribution of such modes, each of which 
is itself a solution that oscillates with some real frequency. We described how Hamiltonian 
systems with continuous spectra generally have a solution formula in terms of such an 
integral over singular eigenmodes. This type of formal expansion led to a an explicit integral 
transformation that transforms the original Caldeira-Leggett system into a pure advection 
problem, just as is the case for the Vlasov-Poisson system. It was noted that a general 
class of such transformations was written down in ji[ and was subsequently extended to a 



larger class of Hamiltonian systems [20| . The existence of these transformations amounts to 
a theory of normal forms for systems with a continuous spectrum, analogous to the theory 
of normal forms for finite degree-of-freedom Hamiltonian systems. This enabled us to write 
down an explicit transformation that converts the time evolution operator for the Caldeira- 
Leggett Hamiltonian into a multiplication operator and we found the inverse of this map. 
In this way we showed that the Caldeira-Leggett model shares the same normal form as the 
Vlasov-Poisson system, along with a number of other Hamiltonian systems that occur in 
different physical contexts. 

One reason for investigating Hamiltonian structure is the existence of universal behavior 
shared by such systems. For example, linear Hamiltonian system with the same normal 
form are equivalent. This suggests some further avenues for research. Here we only treated 
the case where the dispersion relation of the Caldeira-Leggett model has no roots with real 
frequency; i.e. spectrum was purely continuous. When there are roots, the spectrum is no 
longer purely continuous and there are embedded eigenvalues, as is known to be the case 



for the Vlasov-Poisson system 31]. Consequently, one obtains a different normal form, one 
with a discrete component, and this and more complicated normal forms could be explicated. 
We expect that there is a transformation that takes Vlasov-Poisson system with embedded 
modes into the Caldeira-Leggett model with embedded modes. Also, finite degree-of-freedom 
Hamiltonian systems are known to have only certain bifurcations of spectra, for example, as 
governed by Krein's theorem. Since there is a generalization of this theorem for Vlasov-like 



systems 32], one could investigate bifurcations in the context of the Caldeira-Leggett model. 
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Another possibility would be to use the tools developed [27| to do statistical mechanics over 
the continuum bath. Lastly, the integral transform we presented is intimately related to the 
Hilbert transform, which is known to be an important tool in signalprocessing. In the same 
vein the integral transform for the Vlasov-Poisson system of Ref. |3| has been shown to be 
a useful experimental tool [33|, 34 j and one could explore experimental ramifications in the 
context of the Caldeira-Leggett model. 
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